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Abstract 



The running coupling is introduced into the equation for propagation of the pomeron in the 
nucleus via the bootstrap relation. The resulting equation coincides with the one obtained in the 
colour dipole formalism by summing contributions from quark- antiquark loops, with a general 
choice of the regularization scheme. 

(N 
> , 

. 1 Introduction 

O 

. Lately a renewed interest has been shown towards introduction of a running coupling into the 

BFKL dynamics. Summation of contributions from the quark- antiquark loops to the evolution 
of the gluon density has been used to restore the full dependence of the coupling on the running 
scale in the color dipole approach [TJ El [3]. It turned out that the obtained kernel for the linear 
BFKL equation essentially coincides with the one which we found many years ago by imposing the 
bootstrap condition necessary for the fulfilment of unitarity [H[5]. In this paper we draw attention 
to the fact that the bootstrap condition in fact allows to derive also the structure of the triple 
pomeron interaction and thus the form of the non-linear BFKL equation describing propagation 
of the pomeron in the nucleus. This is a simple consequence of the possibility to express the 
basic splitting kernel from 2 to 3 and 4 gluons ('the Bartels vertex', [6j) via the basic BFKL 
interaction. The relation between the splitting vertex and reggeized gluon interaction, together 
with the bootstrap, are the necessary ingredients of the pomeron interaction in the reggeized gluon 
approach. Indeed they allow to present all contributions to the 4-gluon amplitude in the standard 
form of a pomeron splitting into two pomerons. Thus preserving these two relations seems to be 
essential for the construction of pomeron interaction with the running coupling. 

As in our earlier paper, we have to stress from the start that introduction of the running 
coupling into the BFKL formalism cannot be made rigorously and uniquely. The formalism admits 
transverse momenta of any magnitude, including very small ones, at which the concept of the gluon 
and its coupling looses any sense. The introduced running coupling has to be artificially continued 
to small momentum values, where it is completely undetermined. The only information one can 
get from the running of the coupling refers to the region of high momenta. It remains to be seen in 
what measure this information depends on the low momentum behaviour, which cannot be fixed 
theoretically in any reliable manner. It is known however that the non-linear BFKL equation, 
unlike the linear one, is not very sensitive to the infrared region. So hopefully this equation has 
a better chance to produce results weakly dependent on the assumed infrared behaviour of the 
running coupling. 
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The paper is organized as follows. The first section is dedicated to the derivation of the triple- 
pomeron vertex with the running coupling introduced via the bootstrap relation. After recalling 
this method to introduce the running coupling, we go along the same steps as in [HE], where the 
vertex was derived in the limit N c — > oo for the fixed coupling. In Section 3 we construct the full 
amplitude with a single triple pomeron interaction, coupling the vertex with three pomerons. This 
result enables us to build the equation for the pomeron in the large nucleus in Section 4. Finally 
in Section 5 we compare our results with those obtained within the colour dipole approach. 

Having in mind application of the formalism to the large nucleus as a target, we limit ourselves to 
the case of propagation of pomerons with zero total momenta. However generalization to non-zero 
total momenta is straightforward. 



2 Triple-Pomeron vertex 
2.1 Generalities 

As mentioned in the Introduction, this paper follows the idea to introduce a running coupling via 
the bootstrap [U [5]. Derivation of the triple-pomeron vertex in the limit N c — > oo then goes as 
presented in |8] for the fixed coupling case. 

Basic formulas for the introduction of a running coupling via the bootstrap condition consist 
in expressing both the gluon trajectory u and intergluon interaction in the vacuum channel K in 
terms of a single function 77(g) of the gluon momentum, which then can be chosen to conform to 
the high-momentum behaviour of the gluon density with a running coupling. Explicitly 

In these definitions it is assumed that q\ + q 2 = q[ + q' 2 = q- For arbitrary rj(q) the following 
bootstrap relation is satisfied: 

1 f d 2 q' 

a K(q,qi,q[) =uj(q) -uj(qi) -cv{q 2 ). (3) 



2 J (2tt) 2 

The fixed coupling corresponds to the choice 

V(Q) = %1 2 - (4) 
T 

Then one finds the standard expression for the trajectory co(q) and 

K{q,qx,q'x) = ^— -V(q, q x , q' x ), 

where V is the standard BFKL interaction. Note that the extra 2tt in the denominator corresponds 
to the BFKL weight l/(2-7r) 3 in the momentum integration, which we prefer to take standardly as 
V(2vr) 2 . 

From the high-momentum behaviour of the gluon distribution with a running coupling one finds 

71 (q) = l- bq Hn^, q 2 »A 2 , (5) 
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where A is the standard QCD parameter and 

6=^(11JV C -2JV». (6) 

As to the behaviour of 77(g) at small momenta, we shall assume 

7/(0) = 0, (7) 

which guarantees that the gluon trajectory u(q) passes through zero at q = in accordance with 
the gluon properties. The asymptotic ([5]) and condition ([7]) are the only properties of 77(g) which 
follow from the theoretical reasoning. A concrete form of 77(g) interpolating between ([7]) and ([5]) 
may be chosen differently. One hopes that the following physical results will not strongly depend 
on the choice. 

Our old derivation in [7J of the triple-pomeron vertex was actually based on the property ((7|) 
obviously valid for , the bootstrap relation and the relation between the Bartels transition vertex 
-K2-+3 and intergluon BFKL interaction V (Eq. (12) in [7]) 

K 2 ^ 3 (l, 2, 3|f, 3') = V{2, 3|f - 1, 3') - V(12, 3|l'3'). (8) 

Here and in the following we frequently denote gluon momenta just by numbers: 1 = gi, 1' = q[ 
etc. Also we use 12 = q\ + q 2 . All the rest conclusions were obtained from these three relations in 
a purely algebraic manner. 

Our idea is that if we define the transition vertex by a similar relation in terms of the new 
intergluon vertex K, Eq. ([2]), then all the derivation will remain valid also for arbitrary 77(g) 
satisfying (|7j) and thus for a running coupling, provided 77(g) is chosen appropriately. In the next 
subsection we briefly recapitulate successive stages of the derivation for arbitrary 77(g) with 7/(0) = 0. 



2.2 Leading order in N c — > 00 

The changes necessary to pass to arbitrary 77(g) are minimal. Obviously in our old formulas we 
have to drop the coupling constant g factors, since now g is provided by I/77. We also drop factors 
N c in the interaction because we prefer to include it into the kernel K. With these comments, the 
two-gluon equation becomes 

S 20 D 2 = D 20 + K 12 D 2 , (9) 

where 

Sao =J -1 -w(2), (10) 
j is the angular momentum variable, 

D 20 = N c (f(0,q)-f( gi ,q 2 )) (11) 

and f(qi,q 2 ) is just the quark-antiquark loop with, say, the gluon with momentum qi coupled to 
the quark and the other, with momentum q 2 , coupled to the antiquark. With the running coupling 
D 2 q does not change its form, but both u in (I10p and K in ([9]) are now given by Eqs. ([I]) and ([2]) 
respectively. 

The three-gluon system exists in two colour states, which differ in the ordering of the three 
gluons along the loop, 123 and 213. For each order the equation is 

S 30 D 3 = D 30 + D 2 _3 + \{K 12 + K 23 + K 31 )D 3 , (12) 
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where now 

3 

s 3 o = j-i-E w (i) ( 13 ) 

i=i 



and 



Dg 23) = -D^ = \r-^(D 20 (2) - D 20 (l) - D 20 (3)\. (14) 



The new element is the term D 2 ^ which corresponds to transitions of the initial two-gluon system 
into the final three-gluon system. This transition is accomplished by the 2 — > 3 vertex W which, 
as explained, we define via the new intergluon interaction by a relation similar to (J8|) 

VF(l,2,3|l',3') =#(2,311' - 1,3') -#(12,3|1'3'). (15) 

Explicitly, in terms of rj 

I I'M- V ^ V<k23 ' > V ^ 7 ?( 123 ) 

— W(l, 2, 3|1 , 3 ) - ^ _ r ^ _ - ^ _ ~ r? ( 1 /) 7? (3 _ 3 /) + V (V) V (3') ' (16) 

Note that, as with a fixed coupling constant, 

W(l,2,3|l',3') = W(3,2,l|3',l'). (17) 
We find in full analogy with [7J 

D 2 ^ = y^W(l,2,3|l',3') <8>D 2 (1',3') = J^W(1,2,3), (18) 

where <g> means integration over the intermediate gluon momentum with weight l/(27r) 2 . 

Next step is to show that Eq. (fT2j) for the three-gluon system is solved by the reggeized zero 
term ansatz: 

Z)f 3) = -Df l 3) = J^(d 2 (2) - D 2 (l) - D 2 (3)), (19) 



that is by (TTij) in which the loops are substituted by the full solutions of the two-gluon equation ([9]) . 
The proof is wholly based on the bootstrap and relation (115H and literally repeats the corresponding 
proof in [7J. 

Passing to the 4-gluon system, in the limit N c — > oo we find two configurations differing by 
the order of gluons along the quark-antuquark loop: 1234 and 2134. The equation governing the 
4-gluon system is 

S A0 D 4 = D 40 + D 2 ^ A + D 3 ^ A + -(#12 + # 23 + #34 + #4l)£>4, (20) 

where 

4 

s 4 o = j-i-E w C?)- ( 21 ) 

3=1 



The inhomogeneous terms are 



n (1234) 
^40 , 



±N c (d 20 (1) + D 20 (4) - D 20 (U)), (22) 
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£>4o 134) = \n c (d 20 (2) + D 20 (3) - £> 20 (12) - D 20 (13)) , (23) 

= -^W(l,23,4), =0 (24) 

(the definition of W(l, 2, 3) is given by the 2nd equality in (fTBj) ). 




Repeating the corresponding derivation in [7] we demonstrate that in the limit N c — > oo the 
solution of the 4-gluon equation is again given by the reggeized zero-order terms: 



£,(1234) = ^ C (£> 2 (1) + D 2 (4) - D 2 (U)) (27) 

and 

^(2134) = 1^(^(2) + ^(3) - D 2 (12) - D 2 (13)) . (28) 

The proof is purely algebraic and is wholly based on the bootstrap and relations (fT5|) and ([7]) valid 
for any choice of function 77(g) with r/(0) = 0. 



2.3 The triple pomeron configuration 

Next step is to study the next-to-leading configuration in N c — ► 00 corresponding to the triple 
pomeron interaction. Again the derivation practically literally repeats our old one for a fixed 
coupling constant. The governing 4-gluon equation is similar to Eq. (|20p 

S 40 D 4 = D m + L> 2 ^4 + £>3^4 + £>4^4 + (K 12 + K U )D 4 . (29) 

The formal difference is in the absence of interaction between the two final pomerons, which are 
assumed to be made of gluon pairs (1,2) and (3,4), doubling of the rest interactions acting in the 
vacuum colour channels and in the appearance of the term D^±, which describes transitions from 
the leading colour configuration to the subleading one corresponding to two pomerons. 
The four inhomogeneous terms are 

1 4 4 
^40 = ^(E^o(j) " £Z?20(1J)), (30) 

3=1 3=2 

D 2 _4 = -W(l,23,4), (31) 
= JJ- c ( W (^ 2 > 3|1', 3') ® D? 3 4) (1', 3', 4) - W(l, 2, 4|l', 4') ® Z)f 4) (l', 3, 4') 
+ W(2, 3, 4|2', 4') (1. 2 ', 4) - 3, 4|l', 4') ® L>!; 124) (l', 2, 4)) (32) 
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and 

Z?4^4 = ^(#23 + #14 " #13 " #24)( J D? 234) " D^), (33) 

where the D^s on the right-hand side are given by Eqs. (|27p and (|28p . As with a fixed coupling 
constant, all terms except D40 can be presented as a result of action of a certain operator Z on the 
two gluon state, so that Eq. (|29p can be rewritten as 

S 40 D 4 = D 40 + Z®D 2 + (K 12 + K 34 )D 4 . (34) 

The explicit form of this operator can be expressed in terms of a function 

G(l, 2, 3) = -W(l, 2, 3) - D 2 (l, 23)(u>(2) - w(23)) - D 2 (12, 3)(w(2) - w(12)). (35) 

Then one finds 

^ ® D 2 = ^ (2G(1, 34, 2) + 2G(3, 12, 4) + G(l, 23, 4) + G(l, 24, 3) + G(2, 13, 4) + G(2, 14, 3) 

-G(l, 3, 24) - G(l, 4, 23) - G(2, 3, 14) - G(2, 4, 13) - G(3, 2, 14) 

- G(3, 1, 24) - G(4, 2, 13) - G(4, 1, 23) + G(23, 0, 14) + G(13, 0, 24)) . (36) 

This formula is identical to the old one in [TJ, but with new expressions for both W and uj in terms 
of function r\. 

2.4 The triple pomeron vertex 

Corresponding to the two inhomogeneous terms in Eq. (|34p its solution is split in two terms, 
the double pomeron exchange term, generated by -D40, and the triple pomeron interaction term 
Z®D 2 . However one can simplify the solution transferring the part of the double pomeron exchange 
term leading in the high-energy limit into the triple interaction part ( [Hi [8] ) . This is achieved by 
separating from the total solution the reggeized -D40 term: 

D 4 = D 40 (D 20 ^ D 2 ) + Dl (37) 

The irreducible part D\ proves to be a pure triple pomeron interaction, which satisfies 

D{ = Y <g> D 2 + (K l2 + K U )D{. (38) 

The derivation again uses only the bootstrap, relation (fT5j) between W and K and property ([7]). 
The explicit form for the new triple-pomeron vertex Y turns out to be 

Y ® D 2 = ~ (G(l, 23, 4) + G(l, 24, 3) + G(2, 13, 4) + G(2, 14, 3) + G(12, 0, 34) 

- G(l, 2, 34) - G(2, 1, 34) - G(3, 4, 12) - G(4, 3, 12)) (39) 

For a fixed coupling constant this expression was found long ago in [9] . In our approach it remains 
true also for arbitrary function r](q) and thus for a running coupling introduced by means of this 
function. 



G 



3 Coupling to pomerons 



3.1 Momentum space approach 

In the momentum space coupling the vertex Y to two ougoing pomerons is straightforward. The 
two pomerons are described by a product P(l, 2)P(3, 4) and all one is to do is to integrate this 
product with Y ® D 2 over the gluon momenta 1,2,3 and 4 with 12 and 34 fixed and 1234 = from 
the momentum conservation. 

As with a fixed coupling constant, there are certain properties of the wave function and the 
vertex which simplify the resulting expression. First we expect that the pomeron wave function 
in the coordinate space P{r\,r 2 ) vanishes if the two gluons are located at the same spatial point: 
P(r, r) = 0. This property is well-known for the BFKL pomeron with a fixed coupling constant 
and is related to the behaviour of (|4|) at q = 0. With a running coupling this behaviour does not 
change, so we expect that the coordinate wave function will continue to vanish at r\ = r 2 . As a 
result the last five terms in the sum (|39j) will give no contribution, since they depend only on the 
sum of the momenta in one of the pomerons and so put the two gluons at the same spatial point 
in it. If the mentioned property of the pomeron wave function is violated with the introduction of 
a running coupling, we still can drop the five last terms once we restrict ourselves to the case when 
the two pomerons are taken at zero total momentum: 12=34=0, which is the only case relevant for 
the pomeron propagation through the nucleus. Indeed direct inspection shows that 

G( qi ,q 2 ,q3)=0, if qi = 0, or q 3 = 0. (40) 

Second, due to the symmetry of the pomeron wave function, each of the four remaining terms in 
(j39|) gives the same contribution. So, as for the fixed coupling, the triple pomeron vertex effectively 
reduces to 

Y® D 2 = 2G(1, 23, 4). (41) 

Coupling this vertex to two forward pomerons and using Eq. (|35h we get an expression for 
the triple pomeron interaction amplitude T in terms of function 77 (suppressing the dependence on 
rapidities) 

( - D (1' !')( \ 

I 2[ ' >\ v (l-l')ri{4-4') 77(1 - V)r)(A') V (l')r](A-A')J 

, Inn ,){ V( U ) *7(1) \ ■ 1 D M ^( »7(14) U 

2 21 ' >\ri{V)ri{U-V) ??(l')r/(l - T) > 2 21 ' v?](l')r/(14 — 1') r]{V)r]{A - V)) J 

(42) 

with 12 = 34 = 1'4' = 0. In this formula actually both D 2 and P are the pomeron wave functions 
for the incoming and out going pomerons respectively. Function D 2 is the amputated wave function 
related to P by the relation 

P(l,-1) = ^-D 2 (l,-1). (43) 



Expression (|42h is rather cumbersome. A simpler expression is obtained in the coordinate 
representation, which will be presently derived. 
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3.2 Coordinate space approach 

We present 

(2n) 2 6 2 (q 12 - qi - q 2 )P(qi,q 2 ) = / d 2 r 1 d 2 r 2 P qi2 (r 1 ,r 2 )e i ^ +i ^ , (44) 

where P gi2 (ri,r 2 ) is the coordinate wave function of the pomeron with the total momentum q\ 2 . 
Similarly we present the second pomeron via P g34 (r 3 ,r4). Finally 

G(qi,q 2 + q 3 ,q4) = J d 2 ^ 2 ^ 2 ^^^^-^^^ r 2 , r 3 ). (45) 
Then integration over q\, ...q^ gives 

(2vr) 2 <5 2 (( Zl2 + q u )T = 2 J dr 2 ri d 2 r 2 d 2 r 3 P qi2 {r u r 2 )P 334 (r 3 , r 2 )G(n, r 2 , r 3 ) (46) 

We have 

P qi2 (ri,r 2 ) = e^ r ^P(r l2 ), (47) 
where ri 2 = r\ — r 2 . So (f4*6"J) becomes 

(27T) 2 ( 5 2 ( f? i 2 + q 3A )T = 2 1 d 2 r 1 d 2 r 2 d 2 r3e^[ 9l2ri+ ^ 2+9M ) r2+934r3 ]p(n 2 )P(r 3 2)G'(r 1 , r 2 , r 3 ). (48) 

At this stage we note that if we drop P(ri 2 ) from the integrand we get 

/" dr 2 rid 2 r2d 2 r3e5 i [9 12ri+ fei2+934)r2+g34r 3 ]p( r32 ) G ( riir , 2)r3 - ) 

= y d 2 r 2 c/ 2 r 3 el i K 912+934 ) r ' 2+934r3 ]p(r 32 )G( 9l2 ,r 2 ,r 3 ). (49) 

If qi2 = then this expression vanishes due to property (|40l) . As a result one can substitute in (|48j) 

P(r i2 ) ^ P(ri 2 ) - P(0), P(r 32 ) ^ P(r 32 ) - P(0). (50) 

So whether P(0) is equal to zero or not, for the forward case one can always make it equal to zero 
by substitution ([50]) . So in the following we assume P(0) = 0. 

As follows from the translational invariance, for the overall zero total momentum G(fi, r 2 , r 3 ) = 
Cf(ri 2 , r 32 ). So taking as integration variables r 2 , r\ 2 and r 32 we finally obtain 

T = 2 y d 2 r 12 d 2 r 32 P(r 12 )P(r 32 )G(n,r 2 ,r 3 ). (51) 

Now we have to calculate G(r%, r 2 , r 3 ). Due to P(0) = we may drop all terms containing 
<5 2 ( r i2) and/or 5 2 (r 23 ). The total contribution consists of two parts, the first one coming from the 
term — W in (|35p and the second one from the rest. As for the fixed coupling constant case (see 
[8]) in the first part only the first term in (jl6f) gives a contribution which does not contain <5 2 (ri 2 ) 
nor <5 2 (r 23 ) nor both. Direct calculation gives for this contribution 

Gi(n,r 2 ,r 3 ) = -N c D 2 (n,r 3 ) [ d 2 pfj(p)£(r 12 - p)£(r 32 - p), (52) 
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where fj(r) is the Fourier transform of rj(q) and £(r) is the Fourier transform of 1 /77(g). From the 
four terms in the second part the contribution which does not contain S 2 (r 12) nor <5 2 (r23) nor both 
comes from the first and third terms. Its calculation gives 

G 2 (r 1 ,r 2 ,r 3 ) = ^N c D 2 (r 1 ,r 3 ) J d 2 pf,{p)(f{ri2 - p) + £ 2 (r32 - p))- (53) 

Summing we get the final expression 

G(n,r 2 ,r 3 ) = ~N c D 2 {r 13 )F(n 2 ,r 32 ), (54) 

where 

F(n,r 2 ) = J ' d 2 P f,{p)(i{n- p)-t(r 2 - p)f (55) 

and we also used that for the forward pomeron -D2( r i> r 3) = -D2( r i3)- The triple pomeron contri- 
bution to the amplitude is then obtained as 

T(Y) = N c J dy J d 2 r 12 d 2 r 32 F(r 12 ,r 32 )P(Y - y, r 12 )P(Y - y, r 32 )D 2 (y, r 13 ), (56) 

where we restored the y-dependence, suppressed up to now. 

Note that for a fixed coupling constant we have (see Appendix and also [8]) 

G^r 1 ,r 2 ,r 3 ) = -^^(g 2 D 2 {r 13 )). (57) 

07r '12'32 

We can consider our new expression as a result of changing the fixed g 2 to a running quantity 



47T 3 ? 



2 J2. 



<T-> f^F(r 12 ,r 32 ). (58) 

r 13 

Function F{r\,r 2 ) can be presented in a different form, which demonstrates absence of ultravi- 
olet divergency coming from the singular behavior of fj(p) at p — > 0. We have 

prim in p)-J p {27r)2 {27r)2 (2?r)2 v{qiHq2) e 

d 2 qi d 2 q 2 Viqi + q2)_ eiri ( qi+q2 ) (59) 



(2vr) 2 (2vr)2 rj( qi )rj(q 



2 



and similarly 



2 nf,(nV(r, - nVtr, - n~) = I d 2 n ^ ^ J<lP+m(ri-p)+i<l2(r 2 -p) 



d PVmn - P)t(r2 ~P)= jd p (2?r)2 ^ e 



/ 



d 2 qi d 2 q 2 r l (q 1 +q 2 ) eiriqi+ir2q2 



(2vr)2 (2vr) 2 r,( qi )ri(q 2 
So we find 

^--)-/||||»j(^-« ,5, ™)(^-^)- ^ 

In this form it is clear that F(r\,r 2 ) is a well defined function which does not contain ultraviolet 
nor infrared divergency. 

For further use note the identity 

^^- r , r ^/||^^(^-l)(,---l)/^-"<«.-' =0, (62) 
since r/(0) = 0. 
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4 Pomeron in the nucleus 



With the triple Pomeron vertex known, it is straightforward to obtain the equation which sums fan 
diagrams describing propagation of the pomeron in the nucleus. Repeating the derivation for the 
fixed coupling constant in |10j we find for this sum r) at fixed impact parameter b: 

&(y,r,b) = *i(y,6,r) + -N e / dy' TT d 2 r j 6 2 (r 1 - r 2 + r 3 ) 

2 Jo 

F{r 2 ,r 3 )r ] 2 {-iV 1 )G(y - y', r, n)<£>(y', r 2 , 6)$(y', r 3 , 6). (63) 
Here G(y, r, r') is the pomeron forward Green function satisfying the equation 

(— + H)G(y,r,r')=5(y)r ] ~ 1 (-iV)r ] -\-iV)5 2 (r-r'), (64) 

with the Hamiltonian H for the non-amputated forward wave function given by 

H = 2uj + K\ (65) 
where u is expressed via function r/ according to (pQ) and 

-* % a ff (66> 

The inhomogeneous term <3?i corresponds to the single pomeron exchange: 

$1 (y, r, 6) = ^iV c ^T(6) | d 2 r'G(y, r, r')p(r'), (67) 

where T(b) is the nuclear profile function and p(r) is the color density of the nucleon. Applying 
operator d/dy + H to Eq. (|63p we find the evolution equation for the pomeron in the nucleus 



d 

-+ti)9{y,T,b) = d{y)9 [r,b) + -JS e \ 
where ^ ) o( r ) b), playing the role of the initial condition, is given by 



(^+ J ff)$(y, r, b) = (5(y)* (r, b ) + \ N c J II d 2 r^ 2 (r-r 2 +r 3 )F(r 2 , r 3 )$(y, r 2 , 6)$(y, r 3 , 6), 



$ (r, 6) = ^T(6)r ? - 2 (-iV) / o(r). (69) 

Unlike the case of the fixed coupling constant this equation is not simplified in the momentum 
space. 

To compare with the dipole approach, we rewrite our Eq. (|68p as a whole in the coordinate 
space. To do this we have to rewrite action of the Hamiltonian ([65]) on the wave function in the 
coordinate space. In transforming the pomeron amplitude $(y, q, b) to the coordinate space we 
have to take into account condition 

$(y,r = 0,6) = (70) 

which we have extensively used. Technically it means that we have to add to <&(y,q,b) a term 
proportional to 5 2 {q) which guarantees this property. This is essential to obtain the correct form 
for the linear part of the evolution equation in the coordinate space 



We are greataly indebted to Yu. Kovchegov who pointed out this circumstance. 
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It is convenient to split (|55p into three terms 

F( ri ,r 2 )) = /(n,n) + f(r 2 ,r 2 ) - 2/(n,r 2 ), (71) 

where 

/(n,r 2 ) = | d 2 pfj(p)^ n - p)£(r 2 - p). (72) 
In terms of this function one easily finds 
d 2 c 



(27T) 



2 



e^(g)$(y, g, b) = N c d 2 nf( ri - r, n - r)$(y, n, 6) (73) 



and 

/ (2^ (2X c^W)^ ^ b) = 2N C f cfnfin-r^My^ub). (74) 
Thus in the coordinate space we get 

H$(y, r, b) = -N c J d 2 n (2f(n - r, n) - /(n - r, n - r)) n) + const, (75) 
where const should be taken to ensure property (|70p . As a result we find 

H$(y,r,b) = N c J d 2 r 1 F(r 1 -r,r 1 )$(y,r 1 ,b) = ~N C J d 2 r 1 F(r 1 -r,r 1 )(<S>(y,r 1 ,b)+$(y,r 1 -r,b)^ . 

(76) 

Using identity (1621) we may add to the bracket any function independent of r\, to finally obtain 

H$(y,r,b) = ^N C J d 2 r l F(r l - r, n) ($(y, n, 6) + n - r, b) - $(y, r, 6)) . (77) 

In this form the linear part of the evolution equation avquires the standard colour dipole structure 
(see e.g. [2]) and the whole evolution equation becomes 



d If 
—$(?/, r) = -^ N c J d 2 riF(n -r,n) 



($(y, n,b) + n - r, 6) - $(y, r, 6) - $(y, n, &)$(?/, n - r, 6)) . (78) 
In the limit of the fixed coupling constant we find, dropping the infrared regularization terms, 

/^i,r 2 ) = -^¥§, (79) 



2 2 2' 

where at r\ = r 2 = r one should understand 1/r 2 as regularized in the ultraviolet (see [8]): 

^■JTff+toPWlne, 6-0. (80) 
So in analogy with pQ we may define 3 running coupling constants by 

/(n,r 2 ) = - V' ( ? W ? ) £ S. (81) 
7r z a s (ri,r 2 J rfr| 

with the additional condition a s (r,r) = a s (r), and rewrite Eq. (|78p as 

9 -r / ,s 1 ,, f , 2 ,2 c/ , m ? / N/«s( r i) . a s (r 2 ) o a s (ri)a s (r 2 ) rir 2 

($(y, n, b) + $(y, r 2 , 6) - $(y, r, 6) - n, 6)*(y, r 2 , 6)) . (82) 
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5 Discussion 



Eqs. (|78p and (|82p present our final result for the non-linear BFKL equation with the running 
coupling. We stress that function 77(g) in them is determined only by its asymptotic form (0) 
together with the requirement (|7|). A simple possibility is to choose 

^ q) = 2^ 2ln ( a + jj)' (83) 

with b given by ([6]) and arbitrary a > 1. Also one has to remember that the equations are derived 
only in the leading order in the running coupling. Already subleading terms of the relative order 
1/ ln(g 2 /A 2 ) remain undetermined, since they correspond to the next-to-leading order in the running 
coupling. 

Our final coordinate space equation (|82p fully coincides with Eq. (101) in pQ obtained in the 
dipole formalism (for the forward case). However in our approach the running couplings a s (r) and 
cxsifi, r 2 ) are defined by Eq. (|8ip in a general manner, irrespective to any regularization procedure 
and are determined by the concrete choice of function 77(g). In fact they are fixed only in as far the 
high-momentum behaviour of this function is known and so admit a high degree of arbitrariness. It 
remains to be seen how this arbitrariness influence concrete results which follow from the solution 
of the evolution equation (f82j) . 



6 Appendix. Function F{r 1^2) with a fixed coupling constant 

We check that for a fixed coupling constant function F(ri,r 2 ) indeed passes into the expression 
corresponding to Eq. (|57p . In the fixed coupling case rj(k) is given by Then we find 

2-7T 

V(P) = --V 2 ^) (84) 

and 2 

CM = -^2 (In r-c), (85) 

where c = ln(2/m) + ip{V} and m is the infrared regularizer (gluon mass). Terms containing c and 
thus depending on the infrared regularization cancel in the final result. Performing the integration 
over p with the help of the 5-function we find 



2 

f fix (r 1 ,r 2 ) = -7^3 (On - c)V|lnr 2 + (lnr 2 - c)V 2 2 lnn + 2Vi lnr^ lnr 2 ) . (86) 

V 2 lnr = 27T(5 2 (r), (87) 



(27T) 

Taking into account that 



we find 

f flX = " J^( 2 ^TI? + 2vr < 5 2 (r 1 )(lnr 2 - c) + 2vr5 2 (r 2 )(lnr 1 - c 



1'2 

Forming F(ri, r 2 ) according to (|7ip we see that terms with c cancel and dropping terms proportional 
to 5 2 (ri) or 5 2 (r 2 ) we find 

^(n, r2 ) = -^<^#, (89, 



1'2 



in full correspondence with ()57|) . 
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